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ABSTRACT 

We present the first measurements of the third moment of the voltage fluctuations in a conductor. This tech- 
nique can provide new and complementary information on the electronic transport in conducting systems. The 
measurement was performed on non-superconducting tunnel junctions as a function of voltage bias, for various 
temperatures and bandwidths up to 1GHz. The data demonstrate the significant effect of the electromagnetic 
environment of the sample. 
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1. INTRODUCTION 

Transport studies provide a powerful tool for investigating electronic properties of a conductor. The I(V) charac- 
teristic (or the differential resistance Rdi ff = dV/dl) contains partial information on the mechanisms responsible 
for conduction. A much more complete description of transport in the steady state, and further information on 
the conduction mechanisms, is given by the probability distribution of the current P, which describes both the 
dc current I(V) and the fluctuations. Indeed, even with a fixed voltage V applied, I(t) fluctuates, due to the 
discreteness of the charge carriers, the probabilistic character of scattering and the fluctuations of the population 
of energy levels at finite temperature T. 1 

The current fluctuations are characterized by the moments of the probability distribution P of order two and 
higher. Experimentally, the average over P is obtained by time averaging. Thus, the average current is the dc 
current I = (I(t)), where (.) denotes time average. The second moment (the variance) of P, (SI 2 ), measures the 
amplitude of the current fluctuations, with 5I(t) = I(t) — I. The third moment (SI 3 ) (the skewness) measures 
the asymmetry of the fluctuations *. Gaussian noise P(I) oc exp(— aSI 2 ) is symmetric, so it has no third moment. 
The existence of the third moment is related to the breaking of time reversal symmetry by the dc current; at zero 
bias, 7 = and positive and negative current fluctuations are equivalent, so (SI 3 ) = 0. An intense theoretical 
effort has emerged recently to calculate the third and higher moments of P in various systems. 2-12 However, 
until now only the second moment has been measured in the many systems studied. Interest in the third and 
higher moments has occurred, first, because its character is predicted to differ significantly from that of the 
second moment; in particular, the third moment is insensitive to the sample 's own Johnson noise for the voltage 
bias case, yet is more sensitive, in a very subtle manner, to noise and loading by the environment. Second, 
measurements of the higher moments may provide a new tool for studying conduction physics, complementary 
to the second moment. It is the first set of issues that we address. 

In this article we summarize the first measurements of the third moment of the voltage fluctuations across a 
conductor, (SV 3 ), where SV(t) = V(t) — V represents the voltage fluctuations around the dc voltage V (some of 
the data presented here have been published in Ref. 13). Below we relate this to (SI 3 ). Our experimental setup 
is such that the sample is current biased at dc and low frequency but the electromagnetic environment has an 
impedance ^ 50 SI within the detection bandwidth, 10 MHz to 1.2 GHz. We have investigated tunnel junctions 



*(5I :i ) = ((i" 3 )), the third cumulant of refs. 2-12 



because they are predicted to be the simplest system having asymmetric current fluctuations. However, any kind 
of good conductor can be studied with the techniques we have developed. We studied two different samples, at 
liquid helium, liquid nitrogen and room temperatures. Our results are in agreement with a recent theory that 
considers the strong effect of the electromagnetic environment of the sample. 9 Moreover, we show that certain 
of these environmental effects can be dramatically reduced by signal propagation delays from the sample to the 
amplifier. This can guide future experiments on more exotic samples. 

2. THEORETICAL OVERVIEW 
2.1. PERFECT VOLTAGE BIAS 

We present first the case of voltage bias. The scattering theory applied to transport (the so-called Landauer 
picture) has been very succesful in predicting many properties of phase-coherent systems coupled to reservoirs. 
We thus employ that framework. The two contacts to the sample are reservoirs at temperature T. The sample 
is described by a set of M transmission channels characterized by their transmission probability T n , which we 
take to be energy independent. Except for the quantum point contact, 14 only statistical properties of {T n } can 
be accessed. The case of a diffusive wire is special in that it is well described by random matrix theory, and the 
full probability distribution of the transmission coefficients is known. 15 The tunnel junction of area A with a 
uniform barrier of tranparency t can be modeled by M ~ 2k 2 F A channels, each of which has T n = t <C 1. Until 
now only the first two moments T p = M" 1 ^2 n (T n ) p of the probability distribution of {T n } have been accessed in 
measurements: the conductance G = MT\e 2 /h, which involves the average T\ of the {T„}, and the noise power 
density 5/2 ( m A 2 /Hz) : 

5/2 = er)GVcoth(eV/2k B T) + 2k B T{\ - r^)G (1) 

where 77 = 1 — (T2/T1) is the Fano factor, which involves the variance of {T„}. For an opaque tunnel junction, 
all T„ < 1 and 77 = 1; for a diffusive wire 77 — 1/3. At equilibrium, V = and Eq.(l) reduces to the Johnson 
noise 5/2 = 2k B TG, whereas for eV ^> fc^T (the so called shot noise regime), Eq.(l) leads to 5/2 = et]I t. 

For 5/3 , we extend the case of a voltage biased single channel 4 to obtain for a multichannel conductor: 

5/3 (V, T) = erjGkBT (6t/(W) + (1 - 2t)U) (2) 

(in A 2 /Hz 3 ) with U = eV/k B T and t = (T 2 - %)/{%. - T 2 ). We find that the function f{U) is simply related to 
5/2 by dSji/dV = er/Gf(U). The tunnel junction with uniform barrier of transparency t corresponds to t = t, 
the diffusive wire to t = 2/5. 6,7 We can rewrite Eq. (2) as: 

5/3 (V, T) = &tk B T^ + (1 - 2t) V e 2 GV (3) 

The first term of Eq. (3) is temperature dependent, whereas the second term is simply proportionnal to the 
current and completely temperature independent. 

We consider some limits of Eq. (3). For eV -C k B T ', 5/3 = e 2 i]GV. There is no dependence on temperature. 
For an opaque tunnel junction, 5/3 = e 2 I. This result is also obtained for classical, uncorrelated quantized 
charges of magnitude e passing across a barrier. The correlations due to the Fermion reservoirs are seen only in 
the factor 77, which is less than one only when some T n are finite. The magnitude of 5/3 for eV ~ 2k B T is of 
order erjGk B T, i.e., scales with the temperature. For eV S> k B T, Eq. (3) gives 5/3 = er/G[k B T + eV(l — 2t)]. 
It can be negative when i > 1/2, as for the single quantum channel at large voltage when t > 1/2. Thus, if 
more than half the electrons incident from the left reservoir are transmitted, 5/3 is negative. Still, 5/3 is always 
positive for eV <C k B T, and is positive at all voltages if all T„<1 (in which case 77 ~ 1 and t <C 1). 

It is noticeable that the shape of the 5/2 (V, T) curve does not depend solely on 77, so a careful calibration is 
necessary to obtain 77. In contrast, t can be experimentally extracted from the ratio r of the slopes at high voltage 
and low voltage; r = 1 — 2t. Thus t can be determined even in the absence of calibration of the experiment. 

'''One often finds 5/2 (V = 0) = 4ksTG and S/2(V S> feT/e) = 2e?77 in the literature; the factor 2 corresponds to the 
contributions of positive and negative frequencies. 



The measurement of the third moment of current fluctuations provides information on the distribution of 
the transparencies of the channels of the sample. In particular, since it involves the third power of the T n , it is 
much more sensitive to open channels than the conductance or even than 5/2 . For example, in tunnel junctions, 
the distribution may have a tail that extends close to t — 1 (pinholes). 16 Let us consider a tunnel junction of 
M = 10 8 channels of transparency t = 1CP 5 (i.e., with a resistance about 26f2). This corresponds to T p = t p . If 
one adds to the junction 100 channels of transparency t' = 0.1, the distribution of the channels is modified. The 
average transmission T\, and thus the conductance, is increased by only 1%. The second moment is multiplied 
by 10, but this modifies r\ by only 10~ 3 . In contrast, % is multiplied by 10 6 , so i ~ 10~ 3 , i.e., is multiplied by 
100 ! For a tunnel junction, the temperature dependent part (the first term of Eq. (3) is proportionnal to i and 
thus enhanced by a factor 100, but the second term, proportionnal to (1 — 2t) is not much influenced by the 
additionnal open channels. 

As we show below, our samples correspond to good junction tunnel for which t <~ 10~ 5 . Using the general 
property t < (ry -1 — 1), we deduce from our measurement of r] > 0.95, that t < 0.05. Thus, one can consider that 
our samples arc tunnel junctions with an opaque barrier. For such a sample, the current noise spectral density 
(related below to the second moment) is given by : 5/2 = eGUcoth(eV/2kBT) . Only at high voltage eV » ksT 
does this reduce to the Poisson result 5/2 = el. The third moment of the fluctuations is given by 5/3 = e 2 I. 

2.2. THE NOISE IN FOURIER SPACE 

The quantities Sip are spectral densities, expressed in A p /Hz p_1 . Until now we have supposed that Sip is 
frequency independent, but in general it may depend on frequency. Then the pth moment depends on p — 1 
frequencies /i . . . f p -i- However, it is convenient to express Sip as a function of p frequencies such that the sum 
of all the frequencies is zero. Introducing the Fourier components i(f) of the current, one has, for a classical 
current: 

Sip(h, ...,/ p _i) = . . . i{f P )) 5(h + ■■■ f P ) (4) 

In quantum mechanics, the current operators taken at different times do not commute; they also do not in 
Fourier space, and the question of how the operators have to be ordered is crucial. 3 ' 12 

In the case of the second moment, one has 5/2 (/) = (i(f)i(-f)) = (\i(f)\ 2 )- It measures the power emitted 
by the sample at the frequency / within a bandwidth of 1Hz. This is what a spectrum analyzer measures. 
Experimentally, the current emitted by the sample runs through a series of cables, filters and amplifiers before 
beeing detected (this can be avoided by an on-chip detection 17 ). Thus, the measured quantity is a filtered current 
j(.f) = i{f)d{f) where g{f) describes the filter function of the detection. 

One is often interested in the total power emitted by the sample. This is obtained by measuring the DC 
voltage after squaring j(t), i.e. J j 2 (t)dt. This quantity is related to 5/2 through: 

/ + oo r r + OQ r 

f(t)dt = / / g{h)g{h) WMh)) §(h + h) = / W)\ 2 Si,{j)df (5) 

It is remarkable that the frequency-dependent phase shift introduced by g(f) has no influence, in agreement with 
the fact that 5/2 has the meaning of a power. If the detection bandwidth extend from F x to F 2 , i.e., g(f) = 1 
for Fi < |/| < F 2 and g(f) = otherwise, and if 5/2 is frequency independent between f\ and F 2 , then the 
total noise is given by (SI 2 ) — 25/2 (i<2 — This corresponds to the number of choices of / between F\ and 
F 2 , i.e., F 2 — Fi, multiplied by the number of choices for the sign of /, i.e., 2. 

Let us consider now 5/3 . As we said before, the spectral density of the third moment of current fluctuations 
depend on two frequencies /i and / 2 , and the third moment of the measured current j(t) is given by: 

(SI 3 )= J j 3 (t)dt = j J J g{h)g(f2)g{h)SAhJ2)5(h + f 2 + h)dhdf 2 dh (6) 

We see that now the phase of g(f) matters. More precisely, 5/3 measures how three Fourier components of 
the current can beat together to give a non-zero result, i.e., it measures the phase correlations between these 
three Fourier components. With the same hypothesis as before, for a detection between F\ and F 2 , one has 
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Figure 1. Schematics of the equivalent circuit used for the theoretical model. 

now (SI 3 ) = 3S , /3(F 2 — 2Fi) 2 if F 2 > 2F\ and (SI 3 ) — otherwise. This shows how important it is to make 
a broadband measurement. This unusual dependence of the result on F\ and F 2 comes from the fact that the 
lowest frequency being the sum of two others is 2F\ whereas the maximum frequency one can subtract to that 
in order to have a DC signal is F 2 . The number of choices for /i and f 2 is (F 2 — Fi) 2 /2 and the number of 

possible signs for the frequencies is 6 (all combinations except + + + and ). As we show below, we have 

experimentally confirmed this unusual dependence of (SI 3 ) on F\ and F 2 . 

2.3. EFFECT OF THE ENVIRONMENT 

We now consider the effects of the sample's electromagnetic environment (contacts, leads, amplifier, etc.); the 
sample is no longer voltage biased. The environment emits noise, inducing fluctuations of the voltage across 
the sample, which in turn modify the probability distribution P. Moreover, due to the finite impedance of 
the environment, the noise emitted by the sample itself induces also voltage fluctuations. This self-mixing of 
the noise by the sample is responsible for the Coulomb blockade of high resistance samples, where it induces a 
modification of the I(V) characteristics. 18 We consider the circuit depicted in Fig. 1, at first neglecting time 
delay along the coaxial cable. The noise of the sample of resistance R is modeled by a current generator i. The 
voltage SV is measured across a resistor i?o, which has a current generator io of noise spectral density S i 2. One 
has SV = —Ro(i + io) with Rd = RRo/{R + Ro) (R in parallel with Ro). The spectral density of the second 
moment of the voltage fluctuations is S V 2 = +R 2 D (Si2 + Sq). Thus, the only effect of the environment on the 
second moment of the noise is to rescale the fluctuations (by R 2 D ) and add a bias- independent contribution. 1 In 
contrast, it has been recently predicted that the third moment of P is significantly modified by the environment 9 : 

S v s = ~R 3 D S I3 + 3R%S q ^ + 3R%S I2 ^ (7) 

We now give a simple derivation of this result. One has (SV 3 ) = — i?£,((i 3 ) + 3 (i 2 io) + 3 (Uq) + (if)))- The 
term (ii 2 ,) is zero. The term (i^) is zero for a environment at equilibrium, but can be non-zero if the noise 
that the sample sees comes from an active device like an amplifier. Our detection method is insensitive to this 
term and we will neglect it. The first term on the right of Eq. (7) is like that of the second moment (it comes 
from (i 3 ))- The negative sign results from an increasing sample current giving a reduced voltage. In order to 
understand the other terms, let us first consider the correlator (i 2 io). This correlator measures the effect of the 
external current io on the amount of noise i 2 emitted by the sample. If io varies slowly on the time scale of 
the current fluctuations (i.e., ~ h/(eV)), one can replace i 2 by its time average, Sj2 taken at the instantaneous, 
slowly varying (with a time scale of the order of the RC time of the junction), voltage V(t) = V + SV(t). Since 
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Figure 2. Schematics of the experimental setup. 

SV depends on i , the correlator (i 2 io) docs not vanish. We can calculate it for SV <C fc^T, for which one can 
Taylor expand 5/2 in 5/2 (V + 6V(t)) ~ 5/2 (V) + SV {t)dS p / dV . We deduce (i 2 i ) ~ R D (i 2 ) dSp/dV. This 
leads to the second term on the right of Eq.(7). It represents the effect of the noise emitted by the environment, 
and disappeares if the environment is at zero temperature. Similarly, the correlator (i 3 ) can be seen as the effect 
of the sample's current i acting on the noise i 2 emitted by itself, and leads to the third term on the right of Eq. 
(7) (the time scale are no longer separated, but a rigorous calculation 9 shows that Eq.(7) is nevertheless valid). 
This happens because of the finite impedance of the environment and persists at zero temperature. It represents 
the feedback of the environment. We give below a simple derivation of how to include the effect of propagation 
time in the coaxial cable, which dramatically affects 5y3 . 

Let us remark that there are terms of order SV 2 and higher, that are small and neglected here. Such terms 
affect slightly 5y2 as well. They correspond to photon assisted noise in the low frequency limit Tilu <C fc^T. 1 ' 24 



3. THE EXPERIMENT 
3.1. SIGNAL TO NOISE RATIO 

Let us evaluate the signal-to-noise ratio (SNR) for the measurement of (SI 3 ) as compared to the measurements 
of (I) and (SI 2 ). We will suppose that the noise is dominated by the amplifier's noise current I a, characterized 
by a spectral density 5,4. We note B the bandwidth of the measurement and r the averaging time. We suppose 
the low frequency cutoff to be much smaller than the upper one, such that (SI P ) oc S p_1 . The number of 
statistically independent measurements is N — Br, and the SNR always increases like N 1 / 2 . For the dc current 
the SNR is : 

IN 1 ' 2 IN 



CTl (var/ A ) 1 /2 (2S A B)!/ 2 1 \j 2S A (8) 

where var denotes the statistical variance, var(x) = (x 2 ) — (x) 2 . We suppose the amplifier's noise to be gaussian, 
so var(J^) = a P {(l\)) p = ci p (2SaB) p . a\ = 1, a 2 = 2, a 3 = 15, etc. * As expected, the result for o\ is 
independant of B, since I is a dc quantity. Similarly, one has for the second moment: 

{2eIB)N 1 / 2 

a2 = (varJ 2 )V2 - ^ 0) 

with x 2 = e(eB)/SA- x 2 is the ratio of the shot noise of a current cB (one electron every B~ x seconds) to the 
noise of the amplifier. Finally, the SNR 0-3 of the third moment is: 

(3e 2 /B 2 )iV 1 /2 3 



'The rapid decay of the SNR with p could be avoided by the use of p amplifiers and correlation methods. 



For a sample with high impedance, one should use detection in the audio frequency. Taking r = Is, B = 
10kHz, Sa = 10~ 28 A 2 /Hz, gives x <~ 10~ 3 . This estimate shows how difficult is the measurement of S/3 as 
compared to a conductance, or even a S T 2 measurement. For a sample with an impedance comparable with 
5017, one can use RF technique and have a huge bandwidth. We chose this option. Our setup corresponds to 
B = 1GHz and a noise temperature of ~ 250K (Sa ~ 3. 10~ 22 ), and thus x ~ 3. 10~ 4 . This seems much worse 
than the previous case; however, one can pass much more current through a low resistance sample, which makes 
the two techniques comparable. The use of a cryogenic amplifier, with 10 times smaller noise temperature would 
increase the SNR by 10 3 / 2 - 30. However, one can be limited by the amount of current through the sample, in 
particular if the relevant voltage scale is set by the temperature. 

3.2. THE SAMPLES 

Two samples have been studied. Both are tunnel junctions made of Al/Al oxide/Al, using the double angle 
evaporation technique. 19 In sample A the bottom and top Al films are 50 nm thick. The bottom electrode was 
oxidized for 2 hours in pure O2 at a pressure of 500 mTorr. The junction area is 15 /1m 2 . In sample B, the films 
are 120 nm and 300 nm thick, oxidation was for 10 min, and the junction area is 5.6/xm 2 . 22 

3.3. EXPERIMENTAL SETUP 

We have measured SV(t) 3 in real time (see Fig. 2). The third moment of the voltage fluctuations S V 3 is a very 
small quantity and its measurement requires much more care than Sy2. For example, any non linearity of the 
amplifiers gives rise to a fake signal, and we have checked that our result was independent of the amplitude of 
the signal at the input of the various amplifiers. Moreover, since the signal is very small, much more averaging 
is needed than for S V 2 . 

The sample is dc current biased through a bias tee. The noise emitted by the sample is coupled out to an rf 
amplifier through a capacitor so only the ac part of the current is amplified. The resistance of the sample is close 
to 50 i7, and thus is well matched to the coaxial cable and amplifier. After amplification at room temperature 
the signal is separated into four equal branches, each of which carries a signal proportionnal to SV(t). A mixer 
multiplies two of the branches, giving 6V 2 (t); a second mixer multiplies this result with another branch. The 
length of the cables between the splitter and the mixers has been adjusted to insure that the products correspond 
to the voltage taken at the same time. This has been done with the use of short (Ins long) pulses. The output 
of the second mixer, 8V 3 (t), is then low pass filtered, to give a signal which we refer to as D. Ideally D is simply 
proportional to Sy3, where the constant of proportionality depends on mixer gains and frequency bandwidth. 
The last branch is connected to a square-law crystal detector, which produces a voltage X proportional to the 
the rf power it receives: the noise of the sample (SV 2 ^ plus the noise of the amplifiers. The dc current I through 
the sample is swept slowly. We record D(I) and X(I); these are averaged numerically. This detection scheme 
has the advantage of the large bandwidth it provides (<~ 1 GHz), which is crucial for the measurement. Due to 
the imperfections of the mixers, D contains some contribution of Sy2 , and some contribution of the amplifiers, 
which are independent of I: D = Dq + a^Sy?, + azSyi. When the current is reversed, one observes that X is 
the same, X(—I) = X(I), while D is not, see Fig. 3. This prooves the asymmetry of the current fluctuations. 
We deduce Sy3 = (D(I) — D(—I))/(2as), since Sy2 does not depend on the sign of I § . 

3.4. CALIBRATION 

In order to calibrate the measurement of (<5F 2 ), we measure the Johnson noise emitted by a macroscopic 50f2 
resistor at various temperatures. This gives the gain of the detecion scheme as well as its noise temperature 
(dominated by the first amplifier). Since the third moment of the current fluctuations had never been measured 
before, we cannot rely on any sample. Thus, in order to determine the magnitude and sign of (SV 3 ) we measured 
the signal D when the sample is replaced by a programmable function generator. The output of the generator 
consist of a pseudo-random sequence of voltage steps (at a rate of 10 9 samples per second, lGS/s), the statistics 
of which we calculated to be asymmetric in a known way. By varying the maximum amplitude A of the signal and 
reversing its sign, we can measure the coefficients 012 and 0:3, which vary as A 2 and A 3 respectively. However, 

Hn the case this is not correct, one has to subtract any parasitic contribution, which is possible since the calibration 
procedure gives access to 02 and 03. 
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Figure 3. D vs. X as the voltage across the sample is varied. For opposite voltage, X gives the same result (X( — V) = 
X(V)) whereas D doesn't. This proves the existence of the asymmetry of the fluctuations. 

the actual signal synthesized by the generator is not exactly the one which one specifies. Due to the finite 
output bandwidth of the generator, the steps are rounded, which modifies the statistics of the output signal. We 
measured this statistics with a 20GS/s oscilloscope, and used the result to calibrate our measurement. 

Now that we have checked what SV 3 is f° r our sample, we can use the result to calibrate future experiments. 
In particular we can deduce from the measurement of a tunnel junction what are the parameters characterizing 
the environment (see later), which is essential if we want to study more exotic samples for which the intrinsic 
S13 is unknown. 

4. EXPERIMENTAL RESULTS 

4.1. DIFFERENTIAL RESISTANCE AND SECOND MOMENT OF NOISE 

Sample A was measured at T — 4.2K. Its total resistance (tunnel junction and contacts) is 62.6 £1 The resistance 
of the junction itself is extracted from the fit of S V 2 (V), which is a function of eV/ksT with V the voltage drop 
across the junction. (The small resistance of the contacts adds only a voltage independent contribution to SV 2 )- 
We find Ra = 49.6 f2. We measured the differential resistance of the junction in the range of voltage that 
corresponds to the measurement of SV 3 - Rdiff is voltage independent to within 1%. We extract the Fano factor 
r\ of the sample from the slope of the SV 2 {V) curve at high voltage. One gets r\ = 1 (i.e., the Poisson result, 
corresponding to an opaque junction) with a precision of a few percent. Measurements on similar samples in the 
superconducting state show that the current at finite voltage below the gap A/e follows the BCS prediction, 20 
an indication that the junction has no open channels . 

Sample B was measured at T = 4.2 K, T = 77 K and 290 K. The contact resistance is <~ 1S1 and the resistance 
of the junction is 86 0, independent of temperature. Rdiff shows significant curvature on a voltage scale of 160 
mV (the scale used at 300K and 77K): # di// (-160mV)=83.5 0, # di //(+160mV)=82.5 0. It is almost flat on the 
16mV scale (the scale used at 4 K). This phenomenon is usual in tunnel junction and associated with the finite 
height and asymmetry of the barrier. 21 The Fano factor of this sample is also 77 = 1. Measurements of shot 
noise on similar samples in a different experimental setup lead to the same conclusion. 22 
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Figure 4. Effect of finite bandwidth on the measurement of (5V 3 ). Each point corresponds to a different value of the 
frequencies Fi and F 2 , as indicated in the plot. The data shown here correspond to sample B at T = 77 K. 



4.2. EFFECT OF THE BANDWIDTH ON THE THIRD MOMENT 

As we demonstrated before S T 3 has an unusual dependence on the detection bandwidth, from F\ to F 2 . This 
property is also valid for Sy3, even in the presence of an environment with frequency independent impedance. 
A general expression for that can be found in Ref. 23. 

A powerful check that D really measures SV 3 is given by varying the bandwidth. The scaling of Sy 3 with F\ 
and F 2 (S V 3 oc (F 2 - 2F X ) 2 if F 2 > 2Fi and otherwise) is characteristic of the measurement of a third order 
moment. We do not know any experimental artifact that has such behavior. F\ and F 2 are varied by inserting 
filters before the splitter. As can be seen in Fig. 4, our measurement follows the dependence on {F 2 — 2i 7 i) 2 , 
which cannot be cast into a function of (F 2 — Fi). Each point on the curve of Fig. 4 corresponds to a full 
(£V 3 (/)) measurement (see figures 5 and 6). An asymmetric external noise would just add just a constant to 
D, and thus would not affect our measurement, even though this constant might have the same dependence on 
bandwidth as (SV 3 ). 

4.3. RESULT FOR THE THIRD MOMENT 

Sv^ieV/ksT) for sample A at T = 4.2K is shown in Fig. 5; these data were averaged for 12 days. S V 3 (eV/kBT) 
for sample B at T = 4.2K (top), T = 77K (middle) and T = 290K (bottom) is shown in Fig. 6. The averaging 
time for each trace was 16 hours. These results are clearly different from the voltage bias result (the dash-dotted 
line in Fig. 5). However, all our data are very well fitted by Eq. (7) which takes into account the effect of 
the environment (see the dash lines of Fig. 5 and 6). The unknown parameters are the resistance Rq and 
the environment effective noise temperature Tq (as we will see the latter does not correspond to the real noise 
temperature of the environment) . We checked that the impedance of the samples was frequency independent up 
to 1.2 GHz within 5%. Fig. 5 and 6 show the best fits to the theory, Eq. (7), for all our data. The four curves 
lead to Ro = 42 fi, in agreement with the fact that the electromagnetic environment (amplifier, bias tee, coaxial 
cable, sample holder) was identical for the two samples. We have measured the impedance Z env seen by the 
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Figure 5. Measurement of S V 3(eV/k B T) for sample A (solid line). The dashed line corresponds to the best fit with Eq. 
(7). The dash dotted line corresponds to the perfect bias voltage contribution and the dotted lines to the effect of the 
environment. 

sample. Due to impedance mismatch between the amplifier and the cable, there are standing waves along the 
cable. This causes Z env to be complex with a phase that varies with frequency. We measured that the modulus 
\Z en v\ varies between 30 fi and 70 f2 within the detection bandwidth, in reasonable agreement with R = 42 f2 
extracted from the fits. Ro represents an average over frequency of the impedance seen by the sample. 23 We 
notice that the theory predicts a structure at low voltage for high enough temperature (see Fig. 6, dashed line, 
middle and bottom) , that we observe only at room temperature. We think this comes from the fact that at 77K 
we used a long cable, and thus had a largely frequency dependent impedance Z env , whereas at room temperature 
we used a much shorter cable. The presence of the low voltage structure is indeed very much dependent on the 
value chosen for Ro, and might be washed out by frequency averaging with a long cable. 

We have measured directly the noise emitted by the room temperature amplifier; we find T ~ 100 K. This 
is in disagreement with the parameters T * deduced from the fit of the data. However, we show below that this 
is explained by the finite propagation time along the coaxial cable between the sample and the amplifier. 

4.4. EFFECT ON T * OF THE CABLE PROPAGATION TIME 

To analyze our results, let us consider again the circuit in Fig. 1, a simplified equivalent of our setup. Ro <~ 50 fi 
is the input impedance of the amplifier, which is connected to the sample through a coaxial cable of impedance 
Ro (i.e., matched to the amplifier) . The sample's voltage reflection coefficient is T = (R — Ro)/(R + Ro)- 
In the analysis we present next we neglect the influence of the contact resistance and impedance mismatch of 
the amplifier, but we have included it when computing the theory to compare to the data. The voltage 5V(t) 
measured by the amplifier at time t arises from three contributions: i) the noise emitted by the amplifier at time 
t: R io(t)/2 ; half of i enters the cable, ii) the noise emitted by the sample (at time t — At, where At is the 
propagation delay along the cable) that couples into the cable: (1 — T)Ri(t — At)/ 2 ; iii) the noise emitted by 
the amplifier at time t — 2 At that is reflected by the sample: TRoio(t — 2At)/2 ; thus, 

SV(t) = [*>(*) + r* (t - 2At)} - |(1 - T)i(t - At) (11) 

For At = 0, Eq. (11) reduces to SV = -R D (i + i ) with R D = RR /(R + Ro). Thus, Eq.(ll) for At = 
corresponds to Eq. (7), which is a particular case of Eq. (12b) of Ref. 9. . Thus, (SV 3 ) = -Ri)((i 3 ) +3 (i 2 io) + 
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Figure 7. Measurement of S V 3 (eV/ksT) for sample A at T=4.2 K (solid lines). AO: no ac excitation (same as Fig. 5). 
Al: with an ac excitation at frequency Q/2n such that cos2f2At = +1; A2: cos2f2Ai = — 1; A3: no ac excitation but a 
63 Q resistor in parallel with the sample. The dashed lines corresponds to fits with Eq. (7). 

3 (Uq) + («o)) f° r At = 0. As we already show, the term (i 2 io) leads to the second term on the right of Eq. (7). 
and the term (z 3 ) to the first term of Eq. (7), and, due to the sample noise modulating its own voltage, the 
third term of Eq. (7) as well. 

The finite propagation time does affect the correlator (i 2 io), which has to be replaced by (i 2 (t — At)i (t)). 
Thus the term Sq in Eq. (7) has to be replaced by (TS^ + S io ^ io ^ 2 At)) / '(1 + T), where S io{t ) io{t - 2 At) is the 
spectral density corresponding to the correlator (io(t)io(t — 2 At)}. For long enough At this term vanishes, since 
io emitted at times t and t — 2At are uncorrelated. Thus, the effect of the propagation time is to renormalize 
the noise temperature of the environment T = RoS^/ '(2fce) into T * = T r/(1 + T). 

We now check whether the finite propagation time can explain our values of T *. The cable of length ~ 2m 
corresponds to At being long for the bandwidth we used. As a consequence, the relevant noise temperature to be 
used to explain the data is T *. For sample A, T — 0.11; including the contact resistance and cable attenuation 
one expects T * = 5 K ; for sample B, T = 0.26 and one expects Tq = 21 K. A much shorter cable was used 
for T = 290 K, and the reduction of T is not complete. These numbers are in reasonably good agreement with 
the values of T * deduced from the fits (see Fig. 5 and 6), and certainly agree with the trend seen for the two 
samples. Clearly T * <C T for the long cable. 

If there are some open channels in the tunnel junction, Sj3 is modified according to Eq. (3). Since we know 
that t is small, the effect on the temperature independent part of 5y3 is small. However, the term proportionnal 
to i might influence It appears that this term is indistinguishable from the effect of the noise of the 

environnement (the correlator (i 2 i^). It adds a (negative) contribution to To, given by 5Tq = —tT(l + Rq/R). 
Taking i = 0.05 (the maximum value compatible with our measurement of 77), one has ST <~ — 0.08T. This effect 
is intrinsic to the sample and is not modified by the propagation delay. As a result, it would give an even better 
agreement between the evaluation of T * and the values deduced from the fit. A precise measurement of t would 



require a quieter environment (i.e., the use of a cryogenic amplifier) and/or a better impedance match between 
the sample and the amplifier. 

5. MODIFICATION OF THE ENVIRONMENT 

In order to demonstrate more explicitly the influence of the environment on Sy 3 we have modified the parameters 
To and Ro of the environment and measured the effect on Sy3 . 

To is a measure of the current fluctuations emitted by the environment towards the sample. Its influence on 
Sya is through the correlator (i 2 io)- This correlator docs not require io to be a randomly fluctuating quantity 
in order not to vanish. So we can modify it by adding a signal A sin ftt to io (with Q within the detection 
bandwidth). That way we have been able to modify T * without changing Ro, as shown on Fig. 7. The term 
Si (t)i (t-2At) oscillates like cos2fiAi, and thus one can enhance (curve AI as compared to AO in Fig. 7) or 
decrease T *, and even make it negative (Fig. 7, A2). The curves A0-A2 are all parallel at high voltage, as 
expected, since the impedance of the environment remains unchanged; Rq = 42 f2 is the same for the fit of the 
three curves. The experiment described here, i.e. measuring the current fluctuations in the presence of a dc plus 
an ac voltage, is called photon assisted noise when S V 2 is measured. 24 In this context, a frequency of the order 
of eV/h has to be used, and the mechanism can be understood in terms of photons. In our case, since the third 
moment of the current fluctuations involves an odd number of current operators, the particle-like aspect of the 
ac excitation is less clear and deserves more study. 

Second, by adding a 63 resistor in parallel with the sample (without the ac excitation) we have been able 
to change the resistance of the environment Ro, and thus the high voltage slope of S*y3. The fit of curve A3 gives 
Ro = 24.8 0, in excellent agreement with the expected value of 25.2 fl (63 f2 in parallel with 42 f2). Since this 
makes the reflection coefficient negative (r = —0.22), the presence of the extra resistor also reverses the sign of 
Tq, as expected. This is not due to a negative temperature, but rather to a negative L. 

6. CONCLUSION 

Our data are fully consistent with a third moment of current fluctuations S/3 being independent of T between 4K 
and 300K when the sample is voltage biased, as predicted for a tunnel junction. The effect of the environment, 
through its noise and impedance, is clearly demonstrated. This is of prime importance for designing future 
measurements on samples with unknown third moment. 
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